Optimal phase covariant cloning for qubits and qutrits 
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We consider cloning transformations of equatorial qubits | ip^) — -i=(|0) + e**| 1)) and qutrits 
| ifj^fi) — -4j(| 0} + e % ^\ 1) + e | 2)), with the transformation covariant for rotation of the phases 4> 
and 6. The optimal cloning maps are derived without simplifying assumptions from first principles, 
for any number of input and output qubits, and for a single input qutrit and any number of output 
qutrits. We also compare the cloning maps for global and single particle fidelities, and we show that 
the two criteria lead to different optimal maps. 



PACS numbers: 03.65.-w 03.67.- 



I. INTRODUCTION 



The impossibility of perfectly cloning unknown quan- 
tum states selected from a nonorthogonal set is a typical 
quantum feature Q], and is the basis of the security of 
quantum cryptography 0,0,0] In fact, the potential 
eavesdropper Eve cannot clone the quantum state trans- 
mitted by Alice, recover it from multiple copies, and re- 
transmit it undisturbed to the receiver Bob. Eve, how- 
ever, can try to realise an approximate cloning 0,0,013 
in an optimal way, maximizing the fidelity of the copies 
with the original state, and this is a possible eavesdrop- 
ping strategy. The eavesdropping strategies that are 
known to be optimal so far are actually based on cloning 
attacks Moreover, quantum cloning allows to 

study the sharing of quantum information among several 
parties and it may be applied also to study the security 
of multi-party cryptographic schemes [P4 . 

Generally the values of the fidelity achieved by opti- 
mal cloning transformations depends on the set of al- 
lowed input states. In particular, higher fidelities can be 
achieved for smaller sets of input states, since the more 
information about the input is given, the better the input 
states can be cloned. More precisely, for group covariant 
cloning 0] — where the set of input states is the orbit 
of a given state under the action of a group of unitary 
transformations — the smaller is the group the higher is 
the fidelity averaged over the input states. 

In this paper we will develop a thorough analysis of 
the cloning map that is optimal for equatorial qubits and 
qutrits without any simplifying assumptions, including 
that of group covariance and the requirement that the 
output of the cloning map has support on a symmetric 
tensor product Hilbert space. As we will see in Sect. 
II, these assumptions can be derived from the form of 
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the fidelity that one wants to maximize. More precisely, 
we will derive with no assumption the optimal quantum 
cloning transformation maximizing the fidelity averaged 
uniformly over all states of a Bloch sphere equator 



IVv> = ^[|o) 



J<t>\ 



1) 



(1) 



where {| 0), 1 1)} represent a basis for a qubit and the pa- 
rameter 4> £ [0, 2ir) is the angle between the Bloch vector 
and the a:-axis, with the equator in the xy plane with the 
Bloch sphere. As we will see, such averaged form for the 
fidelity automatically leads to the optimal cloning covari- 
ant under the abelian group U(l) of phase rotations — 
the so-called phase- covariant cloning After the first 
analysis of Ref. Q , where only the upper bounds for the 
fidelity were derived by exploiting a connection between 
optimal phase covariant cloning and phase estimation, in 
Refs. 01 an d O a value for the fidelity that breached 
the bound given in Ref. [jj was found for the 1 — * 3 
cloning, apparently obtained under the same assump- 
tions [lH, Then in Ref. a cloning transformation 
from an arbitrary number of input copies N to an arbi- 
trary number M of output copies was presented, however, 
it was proved to be optimal only for N = 1. In this pa- 
per we will prove that the cloning maps of Ref. 0] are 
generally suboptimal for N > 1 and we will derive the 
optimal ones for any values of N and M. In the deriva- 
tion of the optimal cloning maps we will use the general 
method designed for group-covariant cloning introduced 
in Ref. |l3j, which exploits the correspondence between 
CP-maps and positive operators. We also extended our 
analysis to the case of phase-covariant cloning for qutrits, 
namely for quantum states with dimension d = 3. Here 
the covariance group U(l) x U(l) is still abelian, and 
describes the rotation of two different phases. 

The paper is organized as follows. In Sect. II we 
describe the general theory of optimal phase-covariant 
cloning, giving the definitions of all relevant quantities 
in the qubit case, since for the qutrit case the treatment 
is strictly analogous. The starting point is the maxi- 
mization of a phase-averaged fidelity, which will lead to 
a phase-covariant CP-map with output on the symmet- 
ric Hilbert space of the output copies. Then the the- 
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ory of group covariant cloning of Ref. [13( is shortly re- 
viewed and specialized to the case of phase-covariance. In 
Sect. Ill we derive the optimal phase-covariant cloning 
for qubits for any number of input and output copies, 
giving the fidelities for all cases. In Sect. IV the same 
derivation is given for qutrits with any number of output 
copies, starting from a single input copy. Finally, in Sect. 
V we conclude with a discussion of the results, and with 
some open problems and future perspectives. 



where 



= ^(|0>+e*|l) 



2)) 



(8) 



denotes an equatorial qutrit state, and in place of we 
will consider the two-phase rotation operator U^fi which 
achieves the phase shift U$,g\ Vv,0') = I ' l P^'+4>,9'+e)- 
Upon defining the rotated map C NM as f°U° ws 



W = Uf M C NM {U^ PM )U$« , (9) 



is equivalent to the identity Cnm 
Since the fidelity f[CNM]{4>) is linear versus the cloning 
map CpfM-, the averaged fidelity in Eq. can be also 
written in the form 



NM\ — 



2tt 



^/[4 M ](o) = /[4 M ](o), 



(10) 



II. OPTIMAL PHASE-COVARIANT CLONING f rom Eq. Q; we see that covariance of the map C NM 

c nm for ever y ( t ) - 

A cloning map is a special kind of quantum channel, i. 
e. a trace-preserving completely positive (CP) map. In 
the cloning case, the CP map C goes from input states 
in Ti. to output states in 7i® M , with the output state in- 
variant under the permutations of the M output spaces. 
More generally, if we have N > 1 identical copies avail- 
able, the map goes from an input state p® N on the in- 
put Hilbert space Tim given by the symmetric subspace 
M = {H® N )+ of the tensor product H® N to the out- 
put space M = Ti® M with M > N, and with the out- 
put state permutation invariant. Actually, as we will 
see in the following, the optimal map itself will have 
the output state restricted to the symmetric subspace 
M = {H® M )+, even though, generally, permutation in- 
variance of the state does not imply that the state has 
support in the symmetric subspace. In the following we 
will denote a cloning map from N to M "copies" as Cnm- 

We want to find a cloning map Cnm which minimizes 
the following averaged fidelity 



where clearly f[C NM ](0) = /[Cnm] (</>), and the averaged 
map C NM is obviously defined as 



ptp _ I Y pip 

u iVM — NM ' 



(11) 



NM 



(2) 



f[C NM ]{4>) = Tr[| ^)(^ \® M C NM {\ </v)(^ D] 



for equatorial qubit input states | Vv) defined as 

I^) = 7f(|0>+^|1»- (3) 

We will call the cloning map phase-covariant if it satisfies 
the following covariance relation 

Cnm{U® n P «uf N ) = Uf M C N M(PM)Ui® M , (4) 

where is the unitary phase rotation operator 

U4 = exp - a z )} , (5) 

°~x,y,z denotes the usual Pauli matrices, and p_\f is any 
state in J\f = (7i® N )+. In particular, according to the 
fidelity in Eq. 0, we will consider only input states of 
the form 



pm = I V'oXV'o 



|®iV 



(6) 



The unitary transformation in Eq. J5J gives the phase 
shift U^ltptj,/) — I ip^'+ij,). For qutrits the situation will 
be analogous, with input states of the form 

W = IVo.oK^oI^, (7) 



Since, by definition, the averaged map C NM is phase- 
covariant, Eq. (| 10p> simply means that the cloning map 
minimizing the averaged fidelity @ must be itself covari- 
ant. Therefore, finding the optimal cloning map Cnm 
which minimizes the fidelity @ is equivalent to find the 
optimal phase-covariant map Cnm which minimizes the 
following fidelity 

f NM = f[C NU ](0) = Tr[| ^oX^o \® M C NM (\ V'oXV'o \® N )]- 

(12) 

Moreover, due to orthogonality with the state 
I V'oXV'o |® M , any component of the output state 
Cnm(\ ipo)(ipo \® N ) which is not supported on the 
symmetric subspace (7i® M )+ will give no contribution 
to the fidelity (|12H . Therefore, there will be always an 
optimal cloning map having output on the symmetric 
space (Ti® M )+, and in the following we can restrict our 
attention to such maps only, and take M = (H® M )+. 
We will also consider for comparison the average single 
particle fidelity 

Fnm = ±Tr[(\ ^o)(^o II®"" 1 + 1| ^0)^0 \1® M ~ 2 + 
+ l^- 1 ]^)^ |)|^o)(^o \® N C NM (\ ^o)(^o \® N )}- 

(13) 

As shown in Ref. it is convenient to study co- 

variant CP maps in terms of invariant positive operators 
which are in one-to-one correspondence with CP maps. 
In the present context this means to consider the positive 
operators Rnm defined as 



(14) 
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where 2jv denotes the identity map over Af = (H® N )+, 
and \I)) is the maximally entangled vector on Af ® Af 



N 



\n = £ 



SN,n) 



I SJV,1 



(15) 



n=0 



{\sN,n}} denoting any orthonormal basis for (?{® N ) + , 
that we conveniently choose as follows 



SAT,-, 



C(N,n)~2 J2 P j \ 00...0111...1), (16) 



N-n 



where {Pj N ^} denote the permutation operators of the N 
qubits, and C(N, n) is the binomial coefficient C(N, n) = 
N\/n\(N — n)\. As shown in Ref. 0], one can see that 
Rn m is a positive operator on AA ® Af, which is in one- 
to-one correspondence with the CP-map Cnm, with the 
trace-preserving condition for the map writing in terms 
of the operator Rnm as follows 



Tr M [R 



NM I 



In- 



(17) 



The map Cnm can be recovered from the positive oper- 
ator Rnm as follows 



Cnm{pm) = Tw[(/a4 <S> p\[)Rnm], 



(18) 



where O t denotes the transposed operator of O with re- 
spect to the same orthonormal basis (|15f) chosen for the 
maximally entangled vector in Eq. Ijl6|) . namely one de- 
fines O t = (0 1 *)* where the complex conjugated O* of the 
operator O is defined as the operator having complex 
conjugated matrix of the operator O with respect to the 
same orthonormal basis (|15fl . Notice that for the partic- 
ular state in Eq. JfjJ), one has p^ = p^f, since \tpo)® N 
has all real coefficients on the basis (|15|) . The covariance 
of the CP map Cnm m terms of the operator Rnm 
becomes the invariance relation 



[RNM,U® M ®(U® N y} = 0, (19) 

and in our case we have simply (Uf N )* = U®£ . Then, 
according to the Schur lemmas the positive operator 
Rnm is given by the following direct sum 



Rnm — 



jRv , 



(20) 



where v runs over all incquivalent unitary irreducible rep- 
resentations (UIR) contained in the reducible one Jjf M (g> 

(U® N )*, with all equivalent representations grouped to- 
gether, and with i?„ denoting any positive operator over 
the space of all representations equivalent to is, with the 
overall constraint of the trace-preserving condition Ijl7|l . 

Our purpose is to find the optimal phase-covariant 
cloning map which maximizes the fidelity in Eq. (jHj) . 
which using Eq. (|18f) can be rewritten in terms of the 
positive operator Rnm as follows 



Inm = ^[{\^){^\® M 



\®N- 



R 



NM\ 



(21) 



The derivation for the case of qutrits will be strictly anal- 
ogous to that of qubits. 



III. OPTIMAL CLONING FOR QUBITS 

Since the phase rotation group is abelian, all UIR's of 
the group are unidimensional. The inequivalent represen- 
tations can be conveniently labeled by the nonnegative 
integer is, corresponding to the invariant spaces of vectors 
where the group action is equivalent to multiplication by 
the phase factor exjp(iiscf>). Therefore, in the reduction of 
the representation U® M £§> U®^ , each UIR equivalent to 
the representation v is spanned by a vector of the type 



\M-j-is,j + v)®\N-j,j), 

j = 0, . . . , mia(N, M -v), is = 0, . . . , M - N, 



(22) 



where | N — j, j) denotes a state of N qubits where N — j 
of them are in state | 0) , while the remaining j are in 
state 1 1). 

We will now look for the optimal transformations, 
namely the transformations that maximize the fidelity 
Jnm- As proved above, we can restrict our attention to 
the symmetric subspace and therefore we will consider 
the equivalent representations corresponding to the sym- 
metric states {| SM.j+u) I s N,j),j = 0, ... min(iV, M — is)}^, 
where v labels the inequivalent representations (is — 
0, M — N). In the evaluation of the fidelity we take 
I V>o) = (| 0) + | l))/V2. The fidelity of the map is made 
of contributions of the form 

Tr[(| VoXV-o \® M + N )(\ s Mtj+v )( SMtk+ „ | ® | s Ntj )(s N , k I)] 

\J C(N, j)C(N, k) yJC{M, j + is)C(M, k + is) . 

(23) 

Each block of equivalent representations labeled by is is 
given by the positive operator 

Rv = E r i fe l s M,j+y){sM,k+v | <S> | SN,])(sN.k \ , (24) 

jk 

where the trace preserving condition for the operator 
Rnm leads to 



M-l 
i/=0 



2 = 0, ...N 



(25) 



Since each single contribution to the fidelity (|2*3*|> is posi- 
tive versus j and k, the operators R v that maximize the 
fidelity have positive elements r^ k and the off diagonal 
terms are as large as possible, i.e. r w - k — r' kk . 

Therefore, the operator R v can be written as a (gen- 
erally non normalized) projector R u = \ r u )(r v |, where 
I »V> = Ej r j I s M,j+v) ® I sjvj), and rj = jrf y 

Let us now explicitly construct the cloning map that 
optimizes the fidelity. We will first consider the simple 
case N = 1. Each term i?„ will therefore give the follow- 
ing contribution to /im 



f( M = Tr[(| ^ )(^o \® M ® | VoX^o \)Ru] 

= ^tt (r VC(M, u) + r^C(M, 1 + is))' 



(26) 
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with f 1M = Y,vfiM- 

For odd values of M the largest contribution to the 
fidelity comes from the projector with D = (M — 
l)/2 because in this case both terms J C(M, v) and 
yj C (M, 1 + v) are equal and are maximized simultane- 
ously. Moreover, this contribution is maximized when the 
values of r and 7*1 are maximized, i.e. for r$ 2 = 
r[ M = 1. In this case the optimal map is given by 

RlM = I T , (M-l)/2)(r(M-l)/2 I (27) 

and the fidelity takes the form 

fiM = 2^tC(M, (M - l)/2) . (28) 

For even values of M the optimization procedure is 
more involved because the coefficients y C(M, v) and 
\J C(M, 1 + v) arc different and cannot be maximized si- 
multaneously by a single value of v. In order to derive 
the form of the optimal map let us first notice that the 
same contribution f^ M in Eq. 1(26(1 is also achieved by 
choosing 1/ = M — v — 1 with Tq = r\ l ^ v ^ 1 and 



(29) 



Therefore we can look at contributions due to maps of 
the form 



R'u — kiRv + Rm-v-i)- 



(30) 



By taking into account the relations 1)29(1 , in this case the 
completeness constraint 1(1711 can be written as 



Ew) 2 +E«) 2 = 2 



(31) 



The optimal map is given by the values of v that give the 
maximum contributions to l(26() . namely for V- = M/2—1 
and v + = M/2. Therefore the optimization problem con- 
sists in maximizing the quantity Tq~ A + B, with the 
constraint (r%~ ) 2 + {r[~ ) 2 = 2 and with A = y/C(M, vJ) 
and B = y/C(M, V- + 1). The solution is given by 



= V2- 



B 



= V2 



A 



VA' 2 +B 2 



(32) 



Therefore, the optimal map -Rim for even values of M is 
given by 



Rim = \{\r v _)(r v _ I + |7v + )(rv + |) , 
with Tq~ and given by Eq. i1-'S2[) . and Tq + 



The fidelity takes the form 



1M 



C(M,M/2-l)C(M,M/2) 
' [C(M,M/2-l) 2 +C(M,M/2)] 2 



(33) 



1 1 



(34) 



Consider now the general case — > M. Each contribu- 
tion fpjM ^° ^ ne fidelity takes the form 



In 



NAI 



1 

2 N+M 



min(Af,M-i/) 

^ r*y/C(N,j)C{M,j + v) 



(35) 



The maximum value is achieved for the representation v 
for which both the terms y/C(N,j) and y/C(M, j + v) 
are maximized at the same time. In fact, Jnm 1S a 
convex function of rj defined on the convex domain 12511 , 
and the maximum is achieved on the extremal points 
Tjj = 1 for some v. This also corresponds to maximize 
the r.h.s. of Eq. 1(3 5|l by adding "coherently" all the 
terms in the sum over j for a single value of v. We 
have then to distinguish different cases: for N odd and 
M odd the simultaneous maximization of y // C(N 7 j) and 

\J C(M, j + v) occurs when j = (N- 1)/2 and v — (M — 
N)/2. In this case the optimal cloning map corresponds 
to = 1, and is described by 



RNM = \ r(M-N)/2)(f(M-N)/2 I • 

The fidelity takes the explicit form 



INM — oW + AJ 



N 



y/C(N,j)C(M, (M - N)/2 + J) 

3=0 



(36) 

2 

•(37) 



An analogous argument and the results given in Eqs. 1(36(1 
and 1(37(1 hold also when M and N are both even. 

Consider now the case of even M and odd N, or vi cev- 
ersa. The two terms y/C(N,j) and y/C(M, j + v) are 
maximized at the same time for the two values v± = 
(M — N ± l)/2. In order to derive the optimal map we 
follow an argument analogous to the case N = 1 dis- 
cussed above. Actually, let us first notice that the same 
contribution f^ M in Eq. ((35(1 is also achieved by choos- 
ing 1/ = M - N - v with = r^zf~ v - Therefore, as in 
the case N = 1 we can look at cloning maps of the form 



R'v — h(Ru + Rm-n-v) 



(38) 



By exploiting the relation 
the completeness condition as 



M-N-v 

N-j 



we can write 



EK) 2 + E(^-i) 2 = 2 ' 3 = 0, -N/2. (39) 

As mentioned above, the greatest contributions to ther 
fidelity are given by the blocks with v± = [M — N±l)/2. 
The optimal cloning map will therefore be of the form 



Rnm = \ (R V - + Rv+ ) , 



with the constraints r„ 



M-N-v- 



'N-j 



(40) 



and (rj ) 2 + 



( r jv-j) 2 = 2 , j = 0, ..N/2. The optimization of the fi- 
delity 1(35(1 with the constraints ((39(1 leads to the following 
solutions 



r v- = /o yfC(M,N-j+u^ 

5 y/C(M,j+v-)+C(M,N-j+v-) ' 

N ~ 3 y/C(M,j+v-)+C(M,N-j + u- ) ' 



(41) 



Let us now consider as a quality criterion to optimize 
the cloning map the optimization of the average single 
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particle fidelity Fnm, defined as 
F NM = iTr[(|V>o)<V>0 ll^ 1 - 1 + l|^o><Vo |1® M ~ 2 + 



M 
«(g>M-l 



(42) 



In this case we assume that the operators Rnm are sup- 
ported on the symmetric subspace (7i® M )+. Notice that 
the last requirement is now an assumption because the 
argument after Eq. I|12f> . valid for the global fidelity /nai, 
does not hold for the average single particle fidelity. In 
this case the fidelity of the map is made of contributions 
of the form 

Trtd^^oll^- 1 !^)^!^) 

X (| S M -j- v ){sM-k-v I ® I S N -j)(SN-k I)] 



- 2 «+l 



K J> 3 ' y/C(M,j+v)C(M,j+v+l) 3 + L ' k 



C(N,j)S jjk + ^C(N,j)C(N,j + l) 



■ y/(M-j-v)(j + v + l)6 j+uk 



(43) 



where we have considered k > j. As in the case of the 
global fidelity /nm, let us start from the case N = 1. 

Each term i?„ will therefore give the following contri- 
bution to Fim 



M 



W) 2 + K) 2 + ±rZrW(M-v)(v + l) 



(44) 

For odd values of M the term yj (M — v)(v + 1) is max- 
imized for v — (M — l)/2. The optimal map, as in the 
case of the optimization of the global fidelity, is given by 
Eq. (123 with r[ M ~ 1)/2 = r [ M ~ 1)/2 = 1. The fidelity in 
this case takes the explicit form 



M 



1 



M+l - 
2M , 



(45) 



For odd values of M , we can argue similarly to the case 
of the global fidelity, and therefore we have to maximize 
the quantity (|4*4*j) with the constraint (r^ ) 2 + (r^ - ) 2 — 2. 
In this case the optimal solution corresponds to r/^ — 
= 1. The form of the optimal map is given by 

Rim = X\r v _)(r v _ | + (1 - A) | r v+ ){r v+ \ , (46) 
with < A < 1, and the fidelity takes the form 



Fim = 



1 



y/M(M + 2) \ 

2M J 



(47) 



The above optimal single particle fidelities are the same 
as the ones reported in Ref. 0] , where cloning transfor- 
mations restricted to the symmetric subspace were stud- 
ied and the optimality of the single particle fidelity was 
proved only for N = 1. 



Consider now the general case N — > M. Each contri- 
bution Fpj M to the fidelity takes the form 



EW 1 

r NM — 2 N + i 



min{N,M—u) 



min(N,M-v)-l 

XT E W+iy/C{NJ)C(N,j + l) 

3=0 



y/(M-j-v)(j + V + l) i 

(48) 

The maximum value is achieve d for the representa tion 
v for which both the terms y/ C(N,j)C(N,j + 1) and 
y/ (M — j — v) (j + v + 1) are maximized at the same 
time. We have then to distinguish different cases: for 
iV" odd and M odd this occurs when j = (N — l)/2 and 
v = (M — N)/2. In this case the optimal cloning map 
corresponds to = 1, and is described by 

Rnm = \ r(M-N)/2)( r (M-N)/2 I ■ 
The fidelity in this case takes the explicit form 



(49) 



Fnm = 



N-l 

WW / y 

3=0 



y/C(N,j)C(N,j + l) 



(50) 



x y/[(M + N)/2 - j][(M - N)/2 + j + 1]. 



The results given in Eqs. (|49|l and 15U|) hold also when 
M and N are both even. Notice that these results are 
in agreement with the ones conjectured in Ref. |l4j | for 
generic N and M. 

Consider now the case where N and M have 
different parity, for example M is even and N 
is odd. The two terms y/C(N,j)C(N,j + 1) and 
yf (M — j — v) (j + v + 1) are maximized at the same 
time for the two values v± = (M— iV±l)/2. The optimal 
cloning map will therefore be of the form 

\{Ru- 



R 



NM 



R 



v+, 



(51) 



with the constraints r„- 



. j - r v N _ - and (rj ) 2 + {r v N _ j ) 2 = 
2 , j = 0, ..N/2 — 1. Therefore, we have to optimize the 
fidelity gBl> with v = V- by taking into account the above 
constraints, namely the quantity 



JV-l 



r NM ~2 



Wm E r^r^C(N,j)C(N,j + l) 



Xy/(M-j-V-)(j + V-+l) . 

(52) 

The forms of the coefficients r* cannot be found in gen- 
eral. As an example we explicitly optimize the fidelity 
for N = 2 and odd M. In this case v_ = (M - 3)/2 and 
the form of the coefficients is given by 



-V2 



-1 , 



y/(M-l)(M+3) 

^/(Af-l)(Af+3) + (Af+l) 2 



(53) 



= V 2 -(^o ) 2 
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and the fidelity takes the explicit form 



2M 



VM 2 +2M-1 
V2M 



(54) 



We can see that in general when TV and M have different 
parity the optimal solutions are not in agreement with 
the optimal transformations conjectured in |14|. 

We want to point out that the fidelity Fnm of the 
above optimal cloning transformations in the limit M — > 
oo coincides with the fidelity of optimal state estimation 
for N equatorial qubits [16j . 

Moreover, we want to stress that the cloning transfor- 
mations that optimize the global fidelity coincide with 
the optimal ones for the single particle fidelity only in 
the cases where N and M have the same parity. 



IV. OPTIMAL CLONING FOR QUTRITS 

In this section we will derive the optimal 1 — ► M 
cloning transformations for equatorial qutrit states 



i 



(|o> 



J<l>\ 



1) 



„>0\ 



2)), 



(55) 



covariant under the group of rotations of both phases <f> 
and 9. Again, since the group is abelian, all UIR's of 
the group are unidimensional, and in a way analogous to 
the case of cloning of qubits, when restricting to output 
states supported on the symmetric subspace (H® N )+, 
the equivalent UIR's are spanned by the vectors 

)|2) (56) 

where v\ — 0, . . . , M — 1 and v 2 — 0, . . . , M — v\ — 1 
label the invariant spaces of the UIR's corresponding 
to multiplication by the phase factor exp(ivi(j) + i^fl), 
and | Sk, p ,q) denotes the normalized symmetric state of k 
qutrits with k — p — q qutrits in state | 0), p in state | 1) 
and q in state | 2) (the state | Sk, P , q ) is a superposition of 
k\/(k — q — p)\p\q\ orthogonal states). 

In this case we will have contributions of the following 
type to the fidelity 

Tr[| ^o,o)(^o,o f M+1 {\ sm,» u » 2 )(sm,„ 1+ i,v 2 | ® | 0)(1 1)] 



= p^x Vf(M, Vl ,U2)y/T(M, I/i + 1, V 2 ) 



where we define 



T(M, vx,v 2 ) = ( M _ Vl ™X)\vi\»2\ 



(57) 



(58) 



Since all the above contributions are positive, we can 
apply the same argument as in the case of qubits 
and consider positive operators of the form R VlM2 — 
I r Vl . V2 ){r Ul ^ 2 I: where 



r^ 2 1^,^)10) 



r\ x ^\ s M , vl +i, V2 ) 



l)+r v 2 1 ' U2 \s M , Vl , V2+ i)\2), 



(59) 



and the trace preserving condition for the operator Rim 
leads to 

E W 2 ) 2 = E W) 2 = E (^ 2 ) 2 - 1 . (60) 



where in each sum V\ and v 2 are constrained to give non 
negative entries in the states in Eq. (|59|l . 

Each operator R Ul ,v 2 gives the following contribution 
to the fidelity 

fui,u 2 = 3AT+T {fa'" 2 \/T{M,vi,V2) 

+rf ' V2 y/T(M, W+T^W) + '" 2 y/T(M, i^TT)) * . 

(61) 

The operator R Vl ^ V2 that gives the highest contribution 
to the fidelity is the one where the values of T(M, vi,v 2 ), 
T(M, vi + 1, v 2 ) and T(M, v\, ^2 + 1) are maximized. This 
is easy to establish in the case of M = 3fc + 1, because the 
three above expressions for T are all simultaneously max- 
imized for v\ — v 2 = k. Therefore, the optimal cloning 
map is given by Rk,k with 7*0 = r\ = T2 = 1. The corre- 
sponding fidelity takes the explicit form 



/m — nM-l T (M, 



Af-l M-v 



(62) 



The cases with M = 3k and M = 3fc + 2 are more 
involved because the three values of T that appear in 
Eq. (|61J) cannot be maximized simultaneously. In order 
to find the form of the optimal maps we will follow an 
argument similar to the case of qubits. Notice first that 
the value of the contribution f Vlil>2 to the fidelity does not 
change by performing any permutation of the basis states 
{| 0), I 1), I 2)} for each of the M + 1 qutrits in the opera- 
tor R Vl „ 2 . This means that the three blocks labelled by 
{vi,u 2 ), (y 2l M-v 1 -v 2 ~\) and [M — v x — v 2 — give 
the same contribution to the fidelity. Therefore, the same 
contribution given by the operator R„ uU2 is achieved also 
by the map 

R\M = \(fivi.v 2 + Rv 2 .M-V\-v 2 -\ + Rm—vi — v 2 — X,vx) ; 

with the following identifications 



1/1, 1/2 v-2 ,M— -v\ — vi — 1 M— V\ —u% — \js\ 



'0 



= r, 



Vi_,u 2 _ u 2 ,M-u-i— v 2 —X _ M-Vt—v 2 —\,v-t 

'1 — '2 — '0 

vi,v 2 _ v 2 ,M-v 1 —v 2 —l _ M—v 1 —v 2 -l,v 1 

'2 — '0 — '1 



(64) 



The completeness constraint Ijl7|l along with Eqs. Ij64(l 
lead to 

E w*) 3 + E w 2 ) 2 + E K"") 2 = 3 ■ ( 65 ) 



If we restrict our attention to the family of cloning trans- 
formations described by R Vl ,u 2 we have to fulfill the con- 
straint 



(66) 
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Let us first consider the case M = 3k. From Eq. i|61|) we 
can see that the representation that contributes mostly to 
the fidelity is the one with v\ = v 2 — k, because one of the 
three coefficients T that appear in (|61|l is maximized and 
the other two take the second possible highest value si- 
multaneously. Therefore, we can maximize the fidelity by 
restricting our attention to the block labelled by v\ and 
v 2 . Moreover, since T(M, v\ + 1, v 2 ) and T(M, v\,v 2 + 1) 
have the same value for v\ — v 2 — k, the expression i|61|) 
is invariant under exchange of the coefficients r{ and 



k.k 



Therefore we can set r 1 



k.k 



k.k 



when we look for 
the optimal solution. The optimization of the contribu- 
tion in Eq. i|61|) with V\ = v 2 = k corresponds to maxi- 



mizing the quantity r^' k A + 2r^' h B, with the constraint 
(r k ' k ) 2 + 2(r k ' k ) 2 = 3 and with A = y/T{M, M/3, M/3) 
and B = y / T(M, M/3 + 1, Af/3). The solution corre- 
sponds to 



k,k ; 



k,k 



=V3-2(rf' fc )», 



k,k 



=v / 3 



(67) 



k.k 



k,k 



and all other nonvanishing coefficients given by Eq. (|64|l . 
The corresponding optimal map is then given by 



Ri 



M 



M/3,M/3+Rm/3,M/3 



i-i+R 



M/3-l,M/3) ' 



(68) 



with the nonvanishing coefficients r^' 2 given by Eqs. 
lEU and 

In the remaining case M = 3fc + 2 the optimization 
argument and the final solution are the same as for M — 
3k. Here we maximize the quantity rg' k A + 2r\' k B, with 
A = T(M, (M - 2)/3, (M - 2)/3) and B = T(M, (M - 
2)/3, (Af — 2)/3 + l). The optimal solution can be derived 
analogously to the previous case. 

As in the case of qubits we can derive optimal maps 
for qutrits by maximizing the average single particle fi- 
delity instead of the global one, and by assuming that the 
operator R is supported on the symmetric subspace. As 
in the case of qubits we will see that the optimal maps 
for the average single particle fidelity are not always the 
same as the ones derived above, where the global fidelity 
was maximized. Actually, in this case the contributions 
to the average single particle fidelity F\m are of the form 



±Tr[(|^o,o)(V>o,o ll^" 1 1 ^0,0)^0,0 I) 
x I s M ){s M ,1/1+1,1/2 I®|0)(1|] 



(M-l)! 



' (M-i/i-i/ 2 )!i/i!i/2! 

9 (M-ui-V2-l)\ui\v2\ V Af! 



(A/-i/i-jy 2 -l)K t/ i + 1 ) ! ' y 2! 
M\ 

^y/(M- Vl -v 2 ){ Vl + l) 
\k M {M - 1/1 - v 2 . v\) , 



(69) 



where we define 

AM^^TrKl^oK^oll^- 1 ) 



hy/p(q + l). 



M 



SM,p,q ){SM ,p-l,g+l |J 

(70) 



The arguments leading to the form (|59|l and to the con- 
straints (|65|) hold also in this case. The contributions to 
the fidelity due to the operators R Vl . V2 are given by 



F„ 



= *[(■ 



a iv 

"1 i u 2 „V\ ,l/ 2 
1 



(M 



2r^r v ^K M {M-v 1 -v 2 ,v 1 ) 



fa 



v\ - v 2 ,v 2 ) 

1,^2)] . 



(71) 



As discussed above, the optimal cloning map corresponds 
to optimizing the coefficients for the block R VliU2 
that gives the maximum contribution l|71|) . In the case 
M = 3k + 1, all the three terms Am that appear in 
the expression l|71|l are optimized at the same time for 
v\ = v 2 = k, and therefore the optimal map has the same 
form as the one found by maximizing the global fidelity. 
The fidelity Fim in this case takes the explicit form 



M 



1 (1 + 2M&) 



3M I 



(72) 



Let us now consider the case of M — 3k. By looking 
at Eq. (|71|) we can see that the maximum contribu- 
tion to the fidelity corresponds to v\ = v 2 = k, be- 
cause one of the three coefficients Am is optimized and 
the other two take the second possible highest value si- 
multaneously. Moreover, since Am{M — v\ — v%,v\) and 
Am(M — v\ — v 2 ^v 2 ) have the same value for v\ = v 2 = k, 
the expression (j71(l the optimal solution corresponds to 
r-y — r 2 ' ■ Therefore, the maximum contribution l|71|l 
corresponds to maximizing the quantity 



n k,k k.k * /. a , j « 

2r ' r{ A A + (r{ ) A 



k,k\2 , 



(73) 



with the constraint ('. 



,fe,fc\2 

) ' 



2(r k ' k ) 2 = 3 and with A A = 



^JM{M + 3)/3M and A B = (M 
solution corresponds to 



3)/3M. The optimal 



k.k 



2(r^) 2 , 



k,k 



k,k 



(74) 



The optimal map has the form (|68ll . with the values of 
the coefficients for eack block fixed according to Eqs. 
(El and (CU). 

In the remaining case M = 3k + 2 the optimization 
argument and the final solution are the same as in the 
M = 3k case, with A A = [M + 4) (M + 1) /3M and 
Ab = (M+1)/3M. Notice that the optimal map coincide 
with the case 1 — > 3 derived in . 

We want to point out that the average single particle 
fidelity in the limit M — > 00 coincides with the fidelity of 
optimal double-phase estimation for a qutrit in the state 

(S9 E3- 
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V. DISCUSSION 

In this paper we derived from first principles the op- 
timal quantum cloning transformations that maximize 
the fidelity averaged uniformly over all equatorial qubit 
states. We have seen that such averaged form for the fi- 
delity automatically leads to the optimal phase-covariant 
cloning. We have then derived the optimal N — > M 
cloning transformation using the method 01 designed 
for group-covariant cloning. We have also considered 
phase-covariant cloning for qutrits, and derived the I — > 
M optimal cloning maps. From our analytical results one 
can see that the fidelities are always larger than those ob- 
tained for the universal cloning |8(. Moreover, the fidelity 
for the qutrit cloning is smaller than the corresponding 
one for the qubit. We also found that the form of the 
optimal cloning maps depends on the criterion adopted 
to assess the quality of the transformation. Actually, we 
showed that the maximization of the global fidelity and 
the maximization of the average single particle fidelity in 
general lead to different solutions. 

We want now to emphasize that the general analy- 
sis performed in Sect. II for optimal phase-covariant 
cloning would be exactly the same for any smaller dis- 
crete phase-covariance group, such as for example the 
discrete group Z4 of 7r/2-rotations that is employed in 
the BB84 cryptographic scheme |2]. Moreover, since the 
averaged fidelity is the same as that of the single state 
whose group orbit generates all possible input states, the 



only feature that can depend on the particular group in 
the following analysis is the irreducibility of the repre- 
sentation Uf M <g> {Uf N )*. This is the same for the full 
rotation group U(l) and for its subgroup Z4 for N = 1 
and M < 2, whereas one may expect a slight improve- 
ment of fidelities for larger N and M. 

Finally, we want to stress that the method used in the 
present paper could be easily generalized to any quality 
criteria — also called cost-function — different from the av- 
eraged fidelity in Eq. (J2J and the single particle average 
fidelity l)42[l. As a matter of fact, the averaging of the 
cost function will always lead to a an optimal cloning 
that is covariant, as long as the cost-function is linear in 
the cloning CP-map. However, it will not be necessarily 
true that the optimal cloning map will have output state 
in the symmetric tensor-product Hilbert space. Actually, 
also in the case of the single particle fidelity we found the 
optimal maps starting from the assumption that the out- 
put state is supported on the symmetric subspace. 
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